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ABSTRACT 

We present detailed observations of the bright short-hard gamma-ray burst 
GRB 090510 made with the Gamma-ray Burst Monitor (GBM) and Large Area Tele- 
scope (LAT) on board the Fermi observatory. GRB 090510 is the first burst detected 
by the LAT that shows strong evidence for a deviation from a Band spectral fitting 
function during the prompt emission phase. The time-integrated spectrum is fit by the 
sum of a Band function with -Epcak = 3.9 it 0.3 MeV, which is the highest yet measured, 
and a hard power-law component with photon index —1.62 it 0.03 that dominates the 
emission below ~ 20 keV and above « 100 MeV. The onset of the high-energy spectral 
component appears to be delayed by ~0.1s with respect to the onset of a component 
well fit with a single Band function. A faint GBM pulse and a LAT photon are detected 
0.5 s before the main pulse. During the prompt phase, the LAT detected a photon with 
energy 30.5^2 6 GeV, the highest ever measured from a short GRB. Observation of this 
photon sets a minimum bulk outflow Lorentz factor, P > 1200, using simple 77 opacity 
arguments for this GRB at redshift z = 0.903 and a variability time scale on the order 
of tens of ms for the ~ 100 keV~few MeV flux. Stricter high confidence estimates imply 
r > 1000 and still require that the outflows powering short GRBs are at least as highly 
relativistic as those of long duration GRBs. Implications of the temporal behavior and 
power-law shape of the additional component on synchrotron/synchrotron self-Compton 
(SSC), external-shock synchrotron, and hadronic models are considered. 

Subject headings: gamma rays: bursts — GRB 090510 — radiation mechanisms: non- 
thermal 
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Introduction 



A difficulty in trying to understand gamma-ray bursts is that, at least in terms of the temporal 
structure of their emission, all GRBs differ. When the overall time scales of the emission are 
considered, however, a pattern does emerge. The durations of the ~ lOOkeV-MeV emission from 
GRBs form a bimodal distribution and hence are divided into two classes, namely the short and 
long duration bursts. The short bursts formally have durations < 2s with typical values around 
~ 0.2 s, whereas the long bursts have a distribution that peaks around ~ 30 s with a tail extending 
to several hundreds of seconds (JKouveliotou et al.lll993l ). 



Although the physics of their gamma-ray emission is not well un derstood, these two classes 
of bursts likely origina. te from distinct types of progenitor systems (JLee fc Ramirez-Ruiz 
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Wooslev k. Bloomll2006l ). Long duration bursts are thought to be produced by the core collapse of 
massiv e stars as evidenced by th e direct association of several nearby GRBs {z < 0.3) with SN Ib/c 
events (jWoosley &: Bloomll2006l ). Consistent with this, the afterglow cou nterparts of loiig dur ation 
GRBs tend to lie in star forming regions of low mass, irregular galaxies (JKocevski et al.ll2009l ). By 
contrast, short duration GRBs have been associated with both ea rly and late type host galaxies, 
in proportions that reflect the underlying field galaxy distribution (|Bergeiil2009l ). In the prevailing 
model for short bursts, they are produced in merger events of a compact binary systems composed 
of two neutron stars or a neutron star /black hole pair and so would tend to originate from older 
stellar populations. 

With the launch and successful operation of the Fermi Gamma-ray Space Telescope, a wider 
observational window has been opened through which a greater understanding of GRBs may be 
obtained. The Large Area Telescope aboard Fermi provides significantly greater energy coverage 
(20 MeV to >300 GeV), field-of-view (2.4 sr) and effective area (8000 cm^ at 1 GeV) than its pre- 
decessor EGRET (JAtwood et al.ll2009l ). Owing to its substantially lower deadtime (26 /is vs 100 ms 
for EGRET), the LAT can probe the temporal structure of even the shortest GRBs. In addition, 
the LAT can localize GRBs with sufficiently high precision to enable follow-up observations by 
Swift and ground-based observatories; and at energies > few GeV, the LAT can distinguish GRB 
photons from background with little ambiguity. The Gamma-ray Burst Monitor, the other science 
instrument on Fermi, comprises an array of 12 Nal scintillators and two BGO scintillators and can 
detect gamma-rays from 8 keV to 40 MeV over the full unocculted sky. The combined capabilities 
of the LAT and GBM enable Fermi to measure the spectral parameters of GRBs over seven decades 
in energy. For sufficiently bright bursts, time-resolved spectral analysis is possible over the entire 
energy range. 

The Fermi observations of the short gamma-ray burst, GRB 090510, take full advantage of the 
GBM and LAT capabilities. The LAT emission shows temp oral structure on time scales as short 
as 20 ms. In addition to the usual Band function component (JBand et al.lll993l ). spectral fits reveal 
a hard power-law component emerging in the LAT band 0.1 s after the onset of the main prompt 
emission in the GBM band. Moreover, a ~ 0.2 s delay is observed between the brightening of the 



~ 200 MeV-GeV emission with respect to the strong count increases in the Nal and BGO. These 
behaviors present severe challenges for emission models of GRBs. 

A photon with energy SO.StIg GeV was detected by the LAT 0.829 s after the GBM trigger. 
This event arrived during the prompt phase and is temporally coincident with a sharp feature in 
the GBM and LAT light curves. Given this energy, the tem poral structure of the burst light curve 
and the known distance to the burst, lAbdo et al.l (|2009al ) have used this photon and its arrival 
time to set limits on a possible linear energy dependence of the propagation speed of photons 
due to Lorentz-invariance violation that would require a quantum-gravity mass scale significantly 
above the Planck mass. Similar to several of the long bursts seen by the LAT, GRB 090510 shows 
a high energy extended emission component that is detected by the LAT as late as 200 s after 
the GBM trigger. In the context of GRB outflow models, the properties of this GeV emission 
and the optical and X-ray afterglow observations by Swift place significant constraints on po ssible 



internal and external shock models for the late time emission of this source (IGhirlanda et al. 
De Pasquale et al.lboinl : iGorsi et al.lbood : bao et al.l[2009l : IXumar fc Barniol Duranll2009bl ). 
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In this paper, we report on and discuss the GBM and LAT observations of GRB 090510 during 
the prompt emission phase. In section 2, we give the basic observational details. In section 3, we 
present a timing analysis, including discussion of its designation as a short-hard GRB and a cross- 
correlation analysis between various energy bands to characterize any energy-dependent temporal 
lags. In section 4, we perform spectral analyses of both the time-integrated and the time-resolved 
data; and we demonstrate the significance of the additional power-law component and how the 
spectra evolves over the course of the burst. In section 5, we derive a lower limit on the bulk 
Lorentz factor given the variability time scales and observations of the highest energy photons, 
describe the constraints imposed on leptonic, hadronic and other models of the prompt emission, 
and discuss implications of the detection of the 30.5 GeV photon for models of the extragalactic 
background light (EBL). Finally, we summarize our results in section 6. 



2. Basic Observations 



On 2009 May 10 at 00:22:59.97 UT (hereafter Tq), both the GBM and LAT instruments 
triggered on GRB 090510 (trigger 263607781). During the first second after the trigger, the LAT 
detected 62 "transient class" even ts with energies > 100 MeV and within 10° of the Swift UVOT 
position (JKuin &: HoverstenI l2009l ) and 12 events above 1 GeV. In the first minute post-trigger, 
the LAT detected 191 events above 100 MeV and 30 above 1 GeV. From spectral fits to the time- 
integrated emission over the time range Tq + 0.5 s to Tq -|- 1.0 s (see section [5]), the fluence of the 
burst is (5.03 ± 0.25) x lO'^ergcm-^ in the lOkeV to 30 GeV band. In the 15-150 keV band, the 
fluence is (4.08 ± 0.07) x 10~'^ergcm"2_ 



Other detectio ns of the promp t emission were niade by Swift (JHoversten et al.ll2009ll. AGILE 
(|Longoet al.lbood ). Konus-Wind (|Golenetskii et ahlbood ^. Suzaku WAM dohmori et ahlbood ). 



and INTEGRAL- ACS. The Swift UVOT instrument measured the positi on of the optical afterglow 
counterpart to be R.A.(J2000), Dec(J2000) = 22'^14™12!5, -26°34^59.2^^ (lKuin fc HoverstenI [2OO9I ) . 
Follow-up optical spectroscopy taken 3.5 days la ter by VLT/F0RS2 m easured [Oil] and H/3 emission 
lines at a common redshift of z = 0.903 it 0.001 (JMcBreen et al.ll2010l ). Using standard cosmological 
parameters {[0,^,0,^^^] = [0.73,0.27,0.71]), this corresponds to a luminosity distance of d^ = 
1.80 X 10^^ cm and implies an apparent isotropic energy of Eiso = (1.08 it 0.06) x 10 ^^ erg flOkeV- 



30 Ge V, To + 0.5 s to Tn-H.Os). The op tical afterglow was also detected by NOT (|01ofsson et al. 



20091 ) and GRO ND (lOlivares et al.ll2009l ). In the radio, however, only upper-limits were obtained 
using the VLA (JFrail k Chandralbood l. 



Th e host galaxy of GRB 090510 was identified as a late-type elliptical or early-type spiral 
galaxy (JRau et al.ll2009l ). in contrast to the dwarf irregular, star- forming galaxies that have been 
observed to harbor lon g duration GRBs. This is c onsistent with the diverse types of hosts identified 
with short GRBs fe.g.. iNakailboOTi : bergeilbood l. 



3. Timing Analysis 

In Figure [H we plot the light curves from the different detectors in various energy bands. 
The upper two panels contain the data from the GBM Nal and BGO detectors in energy ranges 
8-260 keV and 0.26-5 MeV, respectively. Both of these light curves show the precursor event at 
To that caused the GBM trigger. The main part of the emission in both of these light curves 
starts approximately at Tq + 0.5 s, though both light curves have a small but significant feature at 
To -|- 0.4 s. The third panel displays all of the LAT events that passed the on-board gamma-filter and 
have at least one reconstructed track, and the fourth panel shows the light curve for the transient 
class selection at energies > 100 MeV. In the bottom panel, the measured photon energies > 1 GeV 
are plotted versus arrival time. The three LAT events in the > 100 MeV light curve with arrival 
times in the interval To + s to Tq -|- 0.2 s are very likely burst photons. Using the count rate in the 
LAT during the 200 s directly preceding Tq as a measure of the background rate, the probability 
that these three photons would arise by accident is 1.2 x 10~^. 

Classification of GRBs as short versus long bursts is based on the T90 or T50 durations. The 
^90 (^50) duration is defined as the time between accumulating 5% and 95% (25% and 75%) of 
the counts associated with the GRB. For the T gn and T50 durations of GRB 090510, we have 



applied the technique described in iKoshut et al.l (j 19961 ) to determine these quantities using data 
from multiple instruments (see Tabled]). The cause of the wide range in T90 durations that we find is 
illustrated in Figure [21 where we have plotted the cumulative counts for two detectors, GBM/NaI6 
and Swift /BAT, integrated over the energy ranges 50-300 keV and 50-350 keV, respectively. Ideally, 
the selection of the 0% and the 100% plateaus, which designate the onset and the end of the burst 
data accumulation, is unambiguous. However, when large background variations exist, the plateau 
selection is not unique as one can see in the upper panel of Figure [2j For the GBM data, the most 
conservative selection provides T90 = 9.0 s, while alternative detector selections give conservative 



lower values, down to 1.5s (see Tabled]) . The Swift/BAT data (lower panel) have a much lower 
background and allow for a reasonably robust setting of the 100% level and yield Tgg = 4.0 s. 
Similar plots have been obtained for other GBM detectors, as well as for the INTEGRAL-SPI 
and the Suzaku-WAM instruments. We find a narrower range of T^q values, compared to the Tgo 
durations for the same dataset. The latter are very sensitive to 5% background variations, while 
the former are more robust. 

A common feature of long duration GRBs is the trend for their hi gher-energy photons (be- 



low ~ 1 MeV) to arrive before the lower-energy ones (JNorris et al.ll2000l ). By contrast, short du 



ration GRBs typically show no evidence for lags between different energy ranges below IMeV 



([Norris &: Bonnelll 120061 ). Therefore, since the Tgo estimates in Tabled] span the formal 2s divide 
between short and long duration bursts, the presence or absence of the "hard-to-soft" evolution 
seen in long bursts can help determine whether GRB 090510 should be classified as a short burst. 

In order to estimate the energy-dependent lags for GRB 090510, we compute the cross- 
correlation function (CCF) between light curves for different energy bands among the various 
detectors. For these data, we define the CCF as a function of the lag r using 

CCl^(rJ= _ (ij 

Here /(tj) and g{ti) represent the content of the light curves we are comparing for bin i at time 
ti = idt and bin size dt; f and g are the mean values of the two light curves evaluated over the 
time interval considered, {ti^,ti^)] and the number of bins in each light curve is A^ = ^2 — ^i + 1- 
At a given value of r, equation ([1]) is simply the linear correlation coefficient known as Pearson's 
r. If A is large, then r approximately follows a normal distribution with zero mean and standard 
deviation l/yiV. In this case, the probability that a value greater than \r\ = |CCF(r)| can arise 
by chance between two uncorrelated light curves is given by 



P(|r|) = l-erf(|r|0V/2), (2) 

where erf is the error function. Since the amplitude of the variations of the CCF are a strong 
function of the bin size, we instead plot the probability P{t) = P(|CCF(r)|) as a function of the 
lag r. As we show in Figure JS] P{t) is largely insensitive to the bin size and yields probabilities 
that we can easily interpret. 

The emission in the LAT band is suppressed during Tq + 0.5 s to Tq + 0.6 s relative to the 
emission in the GBM band, and that suppression can dominate the CCF probability curves and 
appear as a relative lag. Since we also wish to determine if the individual pulses during the main 
part of the burst are correlated (i.e., where most of the "spiky" structure is seen), we consider two 
time intervals for our CCF studies: Tq + 0.5 s to Tq + 1.5 s and Tq + 0.6 s to Tq -|- 0.9 s. The former 
interval includes essentially all of the prompt phase emission, whereas the latter will be more useful 
for evaluating correlations in the spiky structure. 



In the two leftmost plots in Figure [3l we show the CCF probabilities for the Tq + 0.6 s to 
To + 0.9 s interval. The upper plot compares the Nal data with the BGO data. The most significant 
correlation is at zero lag. Since the Nal data cover energies 8-260 keV and the BGO data cover 
260 keV-5 MeV, this result shows that there is no evidence for an energy-dependent lag in the 
main Band function component of the prompt emission. This is consistent with GRB 090510 being 
a short burst. There is also a secondary minimum at r = —0.075 s, and this value for the lag 
corresponds to the rough separation between the pulses in the Nal and first two pulses in the BGO 
light curves and a shift of the BGO light curves by this amount to later times. In the lower left 
plot, no significant correlation at any lag is seen between the Nal light curve and the LAT data 
above 100 MeV; and we have performed Monte Carlo simulations and confirmed that the lack of a 
correlation is not simply due to the relatively low statistics in the LAT data. If indeed there is no 
such intrinsic temporal correlation, this would strongly suggest a different emission region from the 
Band spectral component, and no coupling between the observed emission from these two regions, 
e.g., via the external Compton mechanism: soft photons from the Band component upscattered to 
high energies by relativistic electrons in a different region. In the two rightmost plots of Figure O 
the results for the Tq + 0.5 s to Tq + 1.5 s interval are shown. The strong correlation at zero lag for 
the Nal versus BGO data is confirmed. In the Nal versus LAT data, the observed lag arises from 
the delayed onset of the > 100 MeV emission by ~ 0.2-0.3 s. 

In addition to the lack of time lags for lower energy photo ns, short GRBs ar e also known to 



be substantially harder than long GRBs. For BATSE bursts, iKouveliotou et al.l (|l993l ) found a 
strong anticorrelation between the Tgo duration and the hardness ratio of the burst. As we will 
show in the next section, GRB 090510 has the highest peak energy ever measured for any kind of 
burst and is undoubtedly one of the hardest GRBs seen. Given this result, the range of Tgo and 
Tso durations we find, and the lack of any temporal lags for lower energy photons, we can safely 
consider GRB 090510 to be a short duration GRB. 

Finally, in Table [21 we give the shortest variability timescales, ty, for the To + 0.6s to To + 0.8s 
and To + 0.8 s to To + 0.9 s intervals. These timescales will be used to compute the minimum 
bulk Lorentz factor implied by the observation of the highest energy (> GeV) photons in those 
epochs (Section [5]). The t^ values are given by the full-width at half-maximum of the shortest pulse 
measured in any of the detectors for a given time interval. 



4. Spectral Analysis 

The spectral analysis included data from the most brightly illuminated GBM/Nal detectors, 3, 
6, 7, 8, &: 9, covering an energy range 8keV-l MeV, and from both GBM/BGO detectors, covering 
200 keV-40 MeV. Since we wish to fit spectra on short time scales, we used the Time- Tagged Event 
(TTE) data with the energy overfiow channels removed. The background for each GBM detector 
was found by fitting the data in 30 s time intervals that preceded and followed the prompt emission 
component by ~ 2 s, using a polynomial function, then extrapolating to the times during the 



10 



burst. Custom-made detector response files for the different GBM detectors were created using 
the Swift /WOT location. The LAT photon data were extracted for energies 100 MeV-200 GeV 
and using an energy-dependent 95% PSF acceptance cone centered on the source, where we have 
increased the acceptance cone radius by adding in quadrature the uncertainty in the Swift /\J\OT 
location. Appropriate for short time scale (<C 1 ks) analyses, we used the "transient" class event 
selection and the corresponding instrument response functions known as P6_V3_TRANSIENT. 
The response matrix files were created using the gtrspgen tool from ScienceTools v9rl5p2. The 
background during the burst was computed by averaging the LAT background over several orbits 
of the spacecraft during epochs when it had the same orbital position and pointing. The joint 
spectral fits were performed using the spectral analysis package RMFIT (version 3.1). 

In Table [3l we show the results of joint spectral fits of the GBM and LAT data for the 
time-integrated and time-resolved data, including the definitions of the intervals used for spectral 
analysis. We define four intervals a, b, c, and d corresponding to Tq -|- 0.5 s to Tq + 0.6 s, Tq + 0.6 s 
to To + 0.8 s, To + 0.8 s to Tq + 0.9 s, and Tq + 0.9 s to Tq + 1. s, respectively. The time-integrated 
spectrum from Tq + 0.5 s to Tq + 1.0 s shows a deviation from the standard Band function and can 
be adequately fit with the addition of a power-law spectral component. 



E \ 
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or a Comptonized spectral component. 



"(-)-(Hiv)-(-^)- 

The Band function consists of two power-laws that are smoothly joined near the peak photon 
energy, Speak- For our analyses, we use the Band function in the form 



niE) = ^(T0(^)"exp(-M^) E<E,, 



(a-^l^peak \°-^ 



^ i£"kevtC) «-P(/5-")(Tooliv)' E>E^ 



(5) 



where E^ = {a - /3)Speak/(2 + a) (IBand et al.lll993l '). 



For the time-integrated data, the peak energy of the Band component is -Epcak = 3-9ib0.3 MeV. 
This is the highest peak energy ever measured in a GRB time-integrated spectrum. The addition 
of the power-law component with a photon index of —1.62 it 0.03 significantly improves the fit by 
more than 5 a compared to a single Band function. This is the first short burst for which such a 
hard power-law component has been measured. The power-law component appears to extrapolate 
to energies well below -Epeak and do minates the Band f unction emission below ~ 20 keV, similar to 



the behavior seen in GRB 090902B (jAbdo et al.ll2009bl ). Figure H] shows the counts spectrum of the 



time-integrated data and the Band function + power-law fit. In Figure [5l we plot this composite 
model as a i/Fi, spectrum and also plot the separate contributions for each component. 
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For the time-resolved spectroscopy, we initially considered the data partitioned into 0.1s time 
bins, starting at Tq + 0.5 s. However, for the analyses we present here, we have combined the data 
in the Tq + 0.6 s to Tq + 0.8 s interval into a single bin in order to have sufficient counts at energies 

> 100 MeV to constrain the fit of the LAT data. The parameters from the spectral fits to the 
selected intervals are given in Table [3l The Band component undergoes substantial evolution over 
the course of the prompt phase, starting out relatively soft with -Epeak ~ 3 MeV, evolving to a very 
hard spectrum with -Epeak ~ 5 MeV, accompanied by the appearance of the power-law component at 

> 100 MeV, and then becoming softer again with .Epeak ~ 2 MeV. The extra power-law component 
hints at a similar sort of soft-hard-soft evolution, but these spectral changes do not appear to be 
commensurate with the Band component evolution (see Figure [5}3). 



5. Discussion and Implications 



The emergence of a distinct high-energy spectral component in the prompt-phase spectrum 
of GRB 090510 establishes that a hard emission component in addition to a Band component is 
found in the short hard class of GRBs. Hard power-law com ponents are also found in three long - 
duration GRBs, n amely GRB 090902B (JAbdo et al.ll2009bl ). GRB 090926A JAbdo et ahlboiOal ). 
and GRB 941017 (JGonzalez et al.ll2OO.ll ). In GRB 090510, the LAT data show that the 7-ray fiux 
above 100 MeV brightens ~ 0.2 s after the start of the bright phase of GBM emission. This behav- 
i or is similar, but o n a shorter timescale, to the delayed onset in the long-duration GRB 080916C 
(JAbdo et al.ll2009cl ) ar id most other Ferm i GRBs, including GRB 081024B, the first short GRB ob- 
served with the LAT (jAbdo et ahlboiObl ). Furthermore, GR B 090510 displays > 100 M eV emission 
significantly extended beyond the duration of the GBM fiux (JDe Pasquale et al.ll20ld ). as observed 
in other Fer mi GRBs and earlie r from GRB 940217 with the EGRET instrument on the Compton 
Observatory (JHurley et al.lll994l ). These behaviors provide important constraints for high-energy 
emission models and could help answer whether the high-energy 7 rays have a leptonic or hadronic 
origin. 

Though the afterglow radiation in b oth long and sh ort GRBs is probably nonthermal syn- 
chrotron emission from an external shock (jSari et al.lll998l ). the situation is less clear in the prompt 
and early afterglow phases when the GRB engine is most power ful. This radiation could be from 



2002- 



Pe'er et al. 



the t hermal photosphere made by the powerful relativistic wind (iMeszaros et al.l 

2007), from magnetic reconnection in Poynting-flux dominated outflows ( Lyutikov fc Blandford 



2003il. or from nonthernial leptonic emissions fo rmed by internal or external shocks (e.g., 
2009bl : bhirlanda et ahlbood : ICorsi et ahlbood ) in the relativistic jet of a GRB. 



Kumar &: Barniol Duran 



De Pasquale et al.l (|20ld ) describe and present models for the Swift and Fermi observations of 



GRB 090510 during the afterglow phase. Here we consider the implications of the prompt phase 
and early afterglow emission for GRB 090510. After deriving the minimum bulk Lorentz factor Fmin 
and considering the various uncertainties that enter into this calculation, we use the observations 
to constrain leptonic synchrotron/SSC model and hadronic models of short duration GRBs. We 
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do not discuss a thermal photospheric interpretation for the Fermi results on GRB 090510. The 
photospheric interpretation overcomes the problem that the GBM spectra are harder than expected 
below -Epeak with the simplest synchrotron emission model (which is only the case at > 2a in 
GRB 090510 during interval b; see Table [3]). Even if it explains much of the GBM emission, 
however, a different origin is needed for the separate hard spectral component observed at LAT 
energies. The coincident narrow spikes between the LAT all events and GBM lightcurves would 
not be easy to explain in a purely photospheric scenario, though Compton-scat tered photospheri c 
emission by internal shocked electrons could produce the coincident components (JToma et al.ll2010l ). 



5.1. Lower Limit on the Bulk Lorentz Factor 



The use of 7-ray observations to constrain the bulk outflow speed of h ighly variable and 



energetic 7-ray emission from GRBs has been stu died by many authors (e.g., iBaring fc: Harding 



1997 



Lithwick fc Sari 



2001 



Razzaque et al.l 12004 ) . A detailed derivation involves an integration 



over the photon spectrum to calculate the opa city of 7 rays emit ted from sources with idealized 
geometries (see Supplementary Information in lAbdo et al.l l2009d ) . A 5- function approximation 
for the 77 opacity constraint gives values of Tram accurate to ~ 10% whenever the target photon 
spectrum is softer than vF^, oc u. In this case, 77 opacity arguments for a 7-ray photon with 
energy ei, in 777,^0^ units , imply a mi nimum bulk Lorentz factor, defined by r^^(ei) = 1, of (e.g. 



Dondi fc Ghisellinilll995l : lNakai]l2007l l. 



min — 



aTdlil + zfhei 



1/6 



2r2 



a + zY^ei 



(6) 



Here fe is the uF^ flux at photon energy nieC^e, which is evaluated at e = e due to the peaking 
of the 77 cross section near threshold. While the local value of the photon index around e has 
some effect on the exact numerical coefficient, this effect is small provided that the target photon 
index is < —1/2. Because of the threshold condition used to relate the high-energy photon and 
the target photons, the solution to equation ^ is iterative but quickly converges. We use this 
expression to estimate T^iin from Fermi observations of GRB 090510 for comparison with more 
accurate calculations. 

For interval b, during which a 3.4 GeV photon was detected, spectral analysis of GBM and 
LAT data during this episode reveals distinct Band- function and power-law components (Figure [5]). 
The Band function has -Epeak = 5.1 MeV, a = —0.48 and /? = —3.09 (Table [3|). The combined 
Band plus power-law fit reaches a peak uFi, flux of ~ 4 x 10~^ erg cm~^ s~^. Writing the variability 
timescale ty{s)= 0.01t_2 s for 10 ms variability timescale, and /^ = 10~^/_5 erg cm~^ s~^, then 
equation 1^ gives Fmin = 1100(^5/t„2)^/® = lO^Fg with ei = 3400/0.511 ^ 6650. The target 
photon energy e = 2r^;^/(l + z)^ei = 110(/-„5/t_2)"'^^^, or ^ 50 MeV, corresponding to the Band 
(3 branch of the function. Depending on whether the 3.4 GeV photon is interacting with the total 
emission or just the photons in the power law, then /_5 ss 0.7 or /_5 ^ 0.1, and T^im — 950 or 
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Fmin = 720, respectively. For interval c from Tq + 0.8 s to Tq + 0.9 s, the same procedure with the 
30.5 GeV photon gives Fmin — 1370 or 1060 for /_5 ^ 0.5 or /_5 ss 0.1 corresponding respectively 
to the combined Band plus power-law fit or the power-law component only. 

The results of numerical integrations to determine Fmin using the more detailed expressions in 



Abdo et al.l (l2009d ) are shown in Table [2j As can be seen, the simple estimates given above are in 
good agreement with the detailed calculation. A number of issues arise in the use of equation ^ or 
the numerical integrations that are important for assessing the value and uncertainty in Fmjn. These 
include the error incurred by the uncertainties in source spectral fitting parameters, which properly 
involves a covariance matrix to correlate uncertainties for different parameters of the Band-function 
and power-law fit. 

For -Emax; we take the highest energy photon associated with that pulse. Table [2] presents the 
values for ty, -Emax, and Fmin for time intervals Tq + 0.6 s to Tq + 0.8 s and Tq + 0.8 s to Tq + 0.9 s, 
which are the only two for which a distinct pulse width could be measured. In interval Tq + 0.6 s 
to Tq + 0.8 s, the Band -|- PL fit shown in Table [3] form the target photon spectrum, while in 
interval Tq + 0.8 s to Tq + 0.9 s, we present results for both the Band and Band -|- PL fits since 
each fits the data reasonably well. We also give values for Fmin assuming that only the hard 
power-law component forms the target photon source. Here we assume that t^ is the same as that 
measured from the BGO emission, which is primarily associated with photons in the Band portion 
of the spectrum. If the variability timescale of the power-law emission is different than assumed, 
which would be compatible with the two components originating from different locations, then 
the minimum Doppler factor limit would change as indicated by equation ([6]). Furthermore, for 
calculations of Fmin we use the spectrum derived on 0.2 s (time interval b) and 0.1 s (time interval 
c) timescales rather than on the shorter variability timescales during which the high-energy photons 
are measured. This is required for accurate spectral analysis, but could underestimate the flux (and 
therefore Fmjn) during the bright narrow spikes, as can be seen from Fig. [TJ 

The derivation of Fmin depends crucially on the assumption that the high-energy radiation 
and target photons are made in the same emitting region. Correlated variability between different 
wavebands would support the cospatial assumption (see Figure [I]), but no strong evidence for 
this behavior was found in the CCF analysis described in Section 3. A conservative assumption 
would be to suppose that the high-energy photon is part of the power-law component and that 
it can potentially interact only with target photons that are part of the same power-law emission 
component. Even in cases where the target MeV photons are made at smaller radii than the high- 
energy photons, or in different regions within the Doppler cone of the emitting surface, spacetime 
overlap will add to opacity, so this should represent the most conservative assumption. Further 
complicating the derivation of Fmin is the assumed emitting geometry and the temporal evolution 
of the radiating plasma. For a blast-wave geometry, the precise value of Fmin depends on whether 
high-energy photons are produced throughout the "shell" or from the inner edge of the "shell," and 



on the dynamical behavior of the target photons (JGranot et al.ll2008l ) 
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Finally, a significant uncertainty on rmin can arise if the photon with observed energy -Emax is a 
random fluctuation of the underlying true spectrum that corresponds to F < Fmin and T^'y(-E'max) ^ 
1. The confidence we have on the value of F depends on the radiative transport and escape of 7 
rays from the emitting region. For interval c from Tq + 0.8 s to Tq + 0.9 s, Fmin = 1218 zt 61 for 
the Band plus power-law fit (Table [2]). Assuming that the intrinsic spectrum extrapolates as a 
power law to high-energies, a likelihood ratio test assuming an exponential escape probability gives 
r/rmin = 0.96, 0.88, and 0.80 and a spherical escape probability gives F/Fmin = 0.89, 0.69, and 
0.49 at the 1, 2, and 3a confidence levels, respectively. The presence of two photons with energies 
between 1 and 2 GeV in interval b and a 7 GeV photon in interval c reduces the likelihood that 
the highest energy photon is a fiuctuation and can be used to independently estimate Fmjn, giving 
a value Fmin ^ 1000. 



GRB 090510 is the second short ORE observed with LAT, after GRB 081024B (lAbdo et al. 



2010bl ). but the first with a redshift, which is required to derive Fmm- The value of Fmin = 1200- 
1300 for GRB 090510 is comparable t o, and slightly lar ger than the values of Fmin — 900 and 
Fmin = 1000 derived for GRB 080916C (|Abdo et al.ll2009cl ^ and GRB 090902B (|Abdo et al.ll2009bl ;i 



using corresponding 77 opacity arguments. This has led to suggestions tha t the GRBs with the most 
luminous LAT emission are those with the largest bulk Lorentz factors (JKumar fc Barniol Duran 
2009bl : iGhirlanda erallbood ). 



5.2. Models for the Prompt Radiation from GRB 090510 

In addition to the requirement of bulk outflow Lorentz factors F > Fmin, models for GRB 
090510 should explain the ~ 0.2 s delay of the onset of the > 100 MeV emission compared to 
the start of the main GBM emission at Tq + 0.5 s, the appearance of a hard component, and the 
high-energy radiation extending to ~To + 150 s. 



5.2.1. Synchrotron/SSC Model 

A standard GRB model for the prompt phase assumes that th e keV-MeV e mission is nonther- 
mal synchrotron radiation from shock-accelerated electrons (e.g., lTavanilll996l ). This emission is 
necessarily accompanied by SSC radiation. 

The SSC component is stronger for a large rat io of nonthermal electron to magnetic-fiel d 
energy density, as expressed by the condition eg 3> es (jSari Sz Esinll200ll : IZhang Sz Meszarosll200ll ). 
and also when the GRB has a lower bulk Lorentz factor for an external shock origin. Here eg 
and cb are the fractions of shocked energy transferred to nonthermal lepton and magnetic-field 
energy, respectively. Lower bulk Lorentz factors that give stronger SSC components result in 
greater attenuation of high-energy 7 rays from 77 pair-production processes, as discussed in Section 
5.1. For generic (-Eiso ~ 10^^ erg, Fq ~ 300, n ~ 1 cm~'^) parameters expected in an external 
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shock model, a hardening of the LAT spectrum due to the deceleration of the blast wave and the 
emergence of the SSC component in the LAT band was expected to take place in the afterglow 



phase, but not in the prompt phase (JDermer et alj 12000 ). With the larger initial Lorentz factors 



Fq > 10^ implied by 77 arguments for LAT GRBs and the earlier emergence of an external shock 
component, an SSC component would persist after the decline of the synchrotron component due 
to less scattering taking place in the Klein-Nishina regime as the blast wave decelerates, and to 
the peak frequency of the SSC flux decreasing from the TeV to the GeV range. The detailed 
observations of GRB 090510 can help determine whether the hard power-law component appearing 
at ~ Tq + 0.7 s can be explained by SSC emission in the LAT waveband during the prompt phase. 

Figure [6] shows results for a numerical model where synchrotron peak energy, peak flux and 
variability time are made to correspond to the observed values shown in FigureO This code employs 
a Compton kernel that accurately treats Compton-scattering of relativistic electrons throughout 
the Thomson and Klein-Nishina regimes, internal 77 opacity, synchrotron self-absorption, radiative 
escape for 7 rays described by exponential and spherical escape probabilities, and second-order SSC. 
The code does not include reprocessing of internally-absorbed radiation or effects of attenuation 
of 7-ray photons by the EBL (see Section 5.3). Parameters appropriate to the Band component 
in interval b (see Figures [1] and [5]) are used, as shown by the dotted curves. The parameters are 
t„ = 14 ms, Epcak = 5.10 MeV, .Emax = 3.4 GeV, a = —0.48, /3 = —3.09, and peak synchrotron flux 
Ri 4 X 10~^ergcm^^ s^^ (see Figure [5] and Tables [3] and [2]) . Figure [6] shows results for F = 500 
(dashed curves) and F = 1000 (solid curves), and for magnetic fields B' = IkG and B' = IMG in 
the upper and lower panels, respectively. Note that subtraction of an underlying hard component 
will not change the peak I'Fi, flux value by more than ~ 5%. 

The unattenuated power of the SSC component is comparable to the synchrotron power when 
B' = 1 kG, whereas the SSC component is much weaker in the strong-field case. The isotropic jet 
power for the B' = I kG model is = 2.0 x lO^'' erg s'^ and = 5.8 x 10^^ erg s"^ for F = 500 and 
F = 1000, respectively, and is dominated by the energy in the escaping radiation. When B' = 1 
MG, the isotropic jet power is = 1.1 x 10^^ erg s"^ and = 7 x 10^'^ erg s"^ for F = 500 and F = 1000, 
respectively, and is dominated by magnetic-field energy. For B' -^ \ kG, the SSC flux becomes 
much brighter than the synchrotron flux, and the jet power becomes dominated by particles, even 
assuming that all particle energy is in the form of relativistic electrons. 

Because the SSC component is strongly attenuated by 77 processes for the F factors considered, 
an electromagnetic cascade will be formed with 7 rays emerging at lower energies where the system 
becomes optically thin. The timescale for the electromagnetic radiation to cascade to energy E^^ is 
shorter than the synchrotron time, given by tgyn ~ QmQ{T /IQmY/'^ /[{B' /'kGfl'^ ^ E^/im MeV] 
s. Unless B' <^ 1 kG, in which case much more power is found in the cascading SSC emission 
than in the synchrotron emission, the cascading timescale is too short to explain the ~ 0.2 s 
delay bet ween the GBM aii d LAT emission. This model also faces the well-known line-of-death 



problem (jPreece et al.lll998l ) that the standard synchrotron mechanism makes a spectrum softer 



than a = —2/3, whereas a = —0.48 it 0.07 in interval b, representing a nearly 3cr discrepancy from 
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the hardest expected synchrotron emissivity. For the strong-field case where the SSC component is 
weak, the separate hard component in GRB 090510 would then have to originate from a different 
mechanism. 



5.2.2. Afterglow Synchrotron Model 



Kumar fc Barniol DuranI (|2009bl lah and Ichirlanda et alJ (|2009l l: Ichisellini eraD (|2009l ) have 



proposed forward shock emission from the early afterglow as the origin of the delayed onset 
and the h ard component ext e nding i nto the LAT ener gy ba nd. This p ossib ility is also con- 



Ghirlanda et al. 



sidered bv iDe Pasauale et al.l (|2010l ). ICorsi et al.l (|2009l ). and iGao et all (|2009l ). In particular 



(J2009l ) calculate the coasting bulk Lorentz factor of the GRB jet by identify- 
ing the time of the peak LAT emission occurring at ~To + 0.7 s with the deceleration time of 
a relativistic blast wave. If r] is the efficiency to convert bulk kinetic energy to 7-ray energy, 
To « 2000 n-V8(ip^^^/o.2 s)-3/8(?7/0.2)-V8(e^ j^^/S.S x 10^^ erg)V8^ ^j^g^.^ n{cm.-^) is the density 
of the surrounding medium. This expression uses an apparent isotropic 7-ray energy release that 
excludes the LAT emission, which if due to synchrotron radiation from relativistic electrons, would 
require a radiative efficiency approaching unity. 

Depending on circumburst density, the implied Lorentz factor is about 2-4 times the value of 
Lmin calculated in the previous section. For this model, the emission radius corresponding to the 
time of the peak LAT flux is i? w 2.4 x lO^^n"^'^ cm. At t w ipeaki the minimum energy electrons 
in the forward shock radiate synchrotron photons at energies 

hv^ ss 3.6(ee/0. l)2^-2(eij/0. 01)^/2 S5^^(t/0. 2 s)-3/2 ]y[gV assuming an electron injection index of 
p ~ 2.4, where ^e is the fraction of the electrons that take part in the non-thermal power-law 
component responsible for the observed emission. 

Recent particle in cell simulations of relativistic collisionless shocks ( Spitkovskvll2008aH bl: lMartins et al 



20091 ) suggest that .^e is fairly small (of the order of a few percent), which may in our case al- 
low hum of several MeV even for low values of es- For eg = 0.1, es = 0.01, and an isotropic 
equivalent ki netic energy of Ek,is o = 10^^ erg (similar to E^^i^o), the cooling break frequency 
is given by (JGranot fc Saril I2OO2I 1. hvc ~ 0.4n-i(ee/0.1)-i(eB/0.01)-i/2^53^/2(t/o.2 s)-V2 keV. 
This expression holds when eb ^ Ce and Klein-Nishina effects are unimportant, so that Y = 
[(1 + Aee/eBY'^ — l]/2 w \J ^el^B > 1- This would imply fast cooling Uc ^ Vm for n ~ 1, which 
could result in a highly radiative shock for eg ~ 1. However, the late time br oad band spectrum 
at t ~ 100 s from optical-UV through X-ray and up to the LAT 7-ray energies ()De Pasquale et al. 



2OIOI ) suggests /iz^c(lOOs) ~ 300 MeV, which, even if with a large uncertainty, together with the 
overall afterglow modeling suggest a much lower external density of n ~ 10~^. Such a low density 
would imply hu^ ^ 0.3{n/10-^)-^{eB/0.0l)-^/'^E~^^'^{t/0.5 s)-^/2QgY (^here the Klein-Nishina 
effect suppresses SSC cooling), so that Vc passes through the (low part of the) LAT energy range 
around ~ 0.5 s from the onset of the main emission episode, or at ~ Tq+1 s, when there is a softening 
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in the LAT photon index, near the end of the prompt emission. Thus, after the prompt emission, the 
LAT energy range would be above both Urn and Uc-, accounting for the observed photon index. This 
would imply huc ~ 18(n/10-5)-i(eij/0.01)-3/2^^3^/2(t/100s)-i/2 j^gV which is consistent with the 
broad band s pectrum at that tirn e, especially since the spectral break around Vc is very smooth 
and gradual (JGranot &: Saril 12002 ). Interestingly, the inferred value of /ifm(lOOs) ~ 0.43 keV and 
the Vm oc t~^/^ scaling gives hum{tpcak ~ 0.2s) ss 4.8 MeV, which is very close to the measured 
value of -Bpeak near tpeak- This model would not produce a spectrum softer than uF^ oc u^'^, so 
would have difficulty accounting for the emission near lOkeV in interval b, which appears to be a 
continuation of the hard spectral component. 

For an adiabatic blast wave and an electron injection index p, the synchrotron flux scales as 
uF^ oc t"~3w/^zy"~P)/2 at frequencies u > Vm^^c- The measured late time LAT flux decay rate of 
uFp oc i-i-38±o.07j^~o.i±o.i ^Q^id [j^ i\^[g QQgQ require p = 2.5 it 0.1, while the spectral slope requires 
p = 2.2 lb 0.2. Both are consistent with p = 2.4 at the la level. This value is also consistent with 



the late time afterglow of GRB 090510 ()De Pasquale et al.ll2010l ). A radiative blast wave at early 
times is not required in order to account for the observed early LAT flux decay rate. 



5.2.3. Hadronic Models 



In hadronic models, photohadronic and proton/ion synchrotron processes induce electromag- 
netic casca des, which lead to synch r otron and Compton emissions from secondar y electron-positron 



Dermer fc Atovanll2006l : iGupta fc Zhangil2007l : lAsano fc Inoud 120071 ). For a target pho- 



oc 



pairs (e.g., 

ton energy distribution n(e) oc e^, the efficiency for photopion processes is oc R~^r~^E~^~ 
T~^t~^E~^~^, where R oc cP^tv is the shock radius. In this expression, protons with ene rgy Ep 
preferentially interact with photons with energies oc V'^/Ep ( WaxmanI 1 1 9951 : iMurasd l2009l ) . The 
large deduced values for F in GRB 090510 make the photopion efficiency low, so that a very 
large energy release is required if the LAT radiation from GRB 090510 is assumed to be from a 
photomeson-induced cascade. 

A stronger magnetic field shortens the acceleration timescale, leading to a larger maximum 
particle energy ii^maxj thus enhancing the photopion production efficiency. In such a str ong magneti c 
field, however, the effective injection index of secondary pairs tends to be about —2 (jCoppilll992l ). 
which yields a flat spectrum in a i^Fi, plot, while the power-law index of the extra component in 
GRB 090510 is ~ —1.6. In a weaker magnetic field, the Compton component from secondary pairs 
can harden the spectrum, though with reduced photopion production efficiency due to the smaller 
value of -Emax- The slower cooling time of protons th an electrons would produce an extended proton- 
induced emis sion feature (IBottcher &: Dermerlll998l ). though the SSC comp onent would d ecay even 



more slowly (IZhang &: Meszaros 



2001 



A recent numerical calculation bv lAsano et al.l (J2009l ) for 
GRB 090510 demonstrates that the proton injection isotropic-equivalent power is required to be 
larger than 10^^ ergs"^ to explain the hard spectra of the extra component in GRB 090510, which 
is ~ 2 orders of magnitude greater than the measured apparent isotropic 7-ray luminosity. 
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An alternative hadron ic scenario is a proton synchrotron model with a very strong magnetic 
field (JRazzaque et al.ll2009l ). This model explains the delayed onset by proton synchrotron radiation 
in the prompt phase due to the time required to accelerate, accumulate, and cool the ultrarelativistic 
protons. For an onset times tonset ~ 0.2 s after the start of the GBM main emission at Tq + O.S s, the 
required magnetic field is B' ^ 7.4 x 10^ {r/1000)-^/''^{tonset/0-'^ s)~2/3(E^/100 MeV)-!/^ q ^^ ^j^g 
shocked fluid frame. The corresponding total energy release, for a two-sided jet beaming factor /^ ~ 
1.5x 10-^(6 Jl den? , is £ ^ 1.3xl05^(r/1000)i6/3(tonsct/0.1 s)5/3(E^/100 MeV)-2/3(0j/l deg)^ erg 



fsee lWanget al.ll2009l . for GRB 080916C). Thus values of T < 1000 and narrow jet angles « 1° 
are required for a proton syn chrotron scen ario, and such strong beaming is not clearly found in 
the short hard class of GRBs (JNakaiJ 120071 ). Proton-dominated GRB models for ultra- high energy 
cosmic rays therefore are plausible only with low values of T and narrow jet collimation to reduce 
the total energy. 



5.3. Implications for the Extragalactic Background Light 

The EBL is dominated by direct starlight in the optical/ultraviolet and by stellar radiation that 
is reprocessed by dust in the infrared. The EBL is difficult to meas ure directly due to contamination 



by zodiacal and Galactic foreground light (JHauser fc: Dweldl200ll ). For sources at sufficiently high 



redshifts, 77 absorption of high-energy 7-rays by EBL photons can provide a means of constraining 
models of the EBL. Figure [7] shows the absorption optical depth, r-y-y, for various models of the 
EBL as a function of 7-ray energy at the redshift z = 0.903 of GRB 090510. 



We have included curves for the two mo dels oflStecker et al.l (20061 ) as well a s the fiducial model 
of Gilmore et al.l (120091 ). the best fit mod el of Kneiske et al.l ( 2004 ). the model bvlFranceschini et al. 
(I2OO8I), and "Model C" of iFinke et al.l (|2009|). The models of Istecker et al.1 (|20od ) border on 
optically thick at 30.5 GeV; all other models considered here give a tr a nsmis sion probability of 
g-T^T, > Q g5_ 'pj^g baseline and fast evolution models of lStecker et al.l (|200a ) give transmission 
probabilities of 0.37 and 0.30, respectively. Although a higher energy ('~ 30.5 GeV) photon was 



found from this burst than for GRB 080916C (13 GeV; lAbdo et al.l (l2009d)), that burst was mor e 
constraining for EBL models due to its higher redshift, z = 4.35 ± 0.15 (JGreiner et al.l l2009l ). 
However, the EBL does evolve with redshift, and this GRB provides an independent constraint 
from a later time and at a closer distance than GRB 080916C. The low optical depth for the 
highest energy photons in GRB 090510 justify neglecting EBL effects in Figure El 



6. Summary 

We have presented Fermi observations of the short hard gamma-ray burst GRB 090510 during 
the prompt emission phase ext ending to 3 s after the G RB trigger (Fermi and Swift observations 
at later times are presented in iDe Pasquale et al.l I2OIOI ) . The apparent isotropic energy measured 
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by the Fermi GBM and LAT from GRB 090510 is (1.08 ±0.06) x 10^^ erg in the energy range from 
10 keV to 30 GeV, where the upper hmit is defined by the highest energy photon. When corrected 
for corresponding energy ranges, this is at least an order-of-ma gnitude grea ter than the isotropic 



energy releases measured from Swift GRBs with known redshift (JNakaril2007l ) , indicating that GRB 
090510 is unusually energetic. A faint precursor and one LAT photon with energy > 100 MeV are 
observed ~ 0.5 s before the main pulse, and two other LAT photons are detected before the start 
of the bright GBM emission. Large numbers of > 100 MeV LAT photons are detected beginning 
at ^To + 0.7s. 

Spectral analysis shows that the time-integrated flux cannot be fit with a single Band spectrum, 
but is well fit with a power-law emission component in addition to the Band component. In the 
time interval Tq -|- 0.6 s to Tq + 0.8 s, the GRB emission can be resolved into a Band function with 
a = —0.48 lb 0.07 and a hard spectral power law with number index = —1.66. The GBM emission 
is observed to vary on a timescale of 14 it 2 ms in this interval, and on a timescale of 12 it 2 ms 
in the time interval c, from Tq + 0.8 s to Tq + 0.9 s. Our results are summarized in Table El In 
time interval c, where the highest energy (~ 30.5 GeV) photon is detected, our most conservative 
assumptions give F > 1000. 

We considered models for the delayed onset of the high-energy LAT radiation and the appear- 
ance of the distinct hard spectral component. In particular, we considered whether the emission 
spectrum can be explained by a synchrotron model for the Band function and SSC emission for 
the hard power-law component in the GeV range. Cascading of SSC emission into the LAT range 
from higher energies could explain the hard component, but the delayed onset is found to be too 
long compared to the time for development of the cascade. Moreover, the hard GBM spectrum 
presents a challenge for synchrotron models that could instead be explained by a photospheric 
emission component. Other models, including external-shock synchrotron and hadronic models, 
are also considered to explain the Fermi observations of GRB 090510. Future Fermi observa- 
tions of short GRBs will help to discriminate between models and improve our understanding of 
relativistic outflows in GRBs. 
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Fig. 1. — In the top four panels, the GBM and LAT hght curves are shown, from lowest to highest 
energies. The bin size for all light curves is 0.01 s. In the bottom panel, the individual photon 
energies as a function of time are plotted. The red vertical dot-dashed line indicates the trigger 
time, and the black vertical dot-dashed lines indicate the boundaries used for the time-resolved 
spectral analysis. These intervals are labeled a, b, c, and d and correspond to Tq + 0.5 s to Tq + 0.6 s. 
To + 0.6 s to To + 0.8 s. To + 0.8 s to To + 0.9 s, and Tq + 0.9 s to Tq + 1. s, respectively. The insets in 
the top four plots show the counts pe r bin within those time intervals. A previous version of this 
figure appeared in lAbdo et al.l (|2009al ). 
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Fig. 2.— GRB 090510 Tgo and T50 estimates using the GBM/Nal 6 (upper panel) and Swift/BAT 
(lower panel) detectors. The horizontal dashed lines indicate the 0% and 100% cumulative plateau 
levels, and the vertical solid and dashed lines indicate the T90 and T50 interval boundaries, respec- 
tively. The arrow in the upper panel indicates the lowest alternative 100% cumulative plateau level, 
which results in T90 = 1.0 s. 
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Fig. 3. — Cross-correlation function probability as a function of the lag r for the specified light 
curve to follow the Nal light curve. The time intervals considered are indicated in the square 
brackets in each panel. The red curves correspond to time bins of size 0.01 s, and the gray curves 
correspond to 0.005 s. The fact that the probability curves are mostly independent of the bin size 
indicates that we are resolving all of the important features of the underlying variability. 
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Fig. 4. — Counts spectrum for the time-integrated (Tq +0.5, To + 1.0 s) data. The Band + power-law 
model has been fit to these data. See Table El 
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Fig. 5. — a. The best-fit Band + power-law model for the time-integrated data plotted as a uFi, 
spectrum. The two components are plotted separately and the sum is plotted as the heavy line. 
The ±1 o" error contours derived from the errors on the fit parameters are also shown, b. The z^F^ 
model spectra (and iblcr error contours) plotted for the different time ranges given in Table O 
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Fig. 6.— Synchrotron/SSC model for interval b (Tq + 0.6 s to Tq + 0.8 s) of GRB 090510. Dotted 
curves show the Band component during this time interval. The magnetic field in the top and 
bottom panels is 1 kG and 1 MG, respectively. The dashed and solid curves show results for T = 500 
and r = 1000, respectively. The upper blue and lower red curves show received fluxes (without 
EEL corrections) assuming spherical and exponential 7-ray escape probabilities, respectively. 
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Fig. 7. — Model predictions of EEL absorption optical depth versus photon energy for GRB 090510. 
The dashed horizontal line indicates an optical depth of unity, and the dashed vertical line shows 
the energy of the highest energy photon. 
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Table 1. Durations of GRB 090510 detected with different instruments. 



Instrument 


Tgo (s) 


T50 (s) 


Energy Range 


GBM/Nal 3 


0.6 


0.2 


50-300 keV 


GBM/Nal 6 


9.0 


0.3 


50-300 keV 


GBM/Nal 7 


1.5 


0.2 


50-300 keV 


GBM/Nal 3, 6 & 7 


2.1 


0.2 


50-300 keV 


Swift/BAT 


4.0 


0.7 


50-350 keV 


INTEGRAL-SPI 


2.5 


0.1 


20keV-10MeV 


Suzaku-WAM 


5.8 


0.5 


50 keV-5 MeV 
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Table 2. Fmin values for the shortest time scale pulses from GRB 090510 





Spectrum 


tv (ms) 


^max (GeV) 


^ mm 


0.6-0.8 


Band + PL 


14 ±2 


3.4 


951 ± 38 


0.6-0.8 


PL^ 


14 ±2 


3.4 


703 ± 34 


0.8-0.9 


Band'^ 


12 ±2 


30.5 


1324 ± 50 


0.8-0.9 


Band + PL 


12 ±2 


30.5 


1218 ±61 


0.8-0.9 


PL^ 


12 ±2 


30.5 


1083 ± 88 



'^Uncertainty on T^nin is associated with uncertainties in ty 
and spectral flux only. 

Variability time tv for power law (PL) component is as- 
sumed to be the same as that derived from the BGO emission, 
which is described primarily by Band component. 

^The Band-function-only fit gives a larger flux than the Band 
plus power-law fit in the portion of the spectrum where the 
highest-energy 7 rays typically interact; see Figure [5l 



Table 3. Prompt emission spectral fit parameters. 



T-T„ 

(s) 



Model 



Band Model 



power-law or Comptonized 



(10- 



A 
-2s-ikcV-i) 



^peak 

(MeV) 



A at 1 GeV 



'pk 



(lO-^cm-^s-ikcV-i) (MeV) 



Index 



CSTAT / dof 



0.5-1.0 


Band 




Band + PL 




Band + Comp 


0.5-0.6 


Band 


0.6-0.8 


Band + PL 


0.8-0.9 


Band 




Band + PL 


0.9-1.0 


PL (LAT only) 



4 S16+°-ll'5 
q 100+0.269 

q r\fiQ-\-0 .281 
0.zuo_Q 266 

2.984l»;fJ 

U.U4U_Q gQ4 

0.028 ± 0.006 



.+0.267 
"'-0.263 

qqR+0.280 
o.yOD_Q 260 

4.UUZ_Q 271 
a+0.185 
"-0.174 
,+0.443 
^-0.400 

1 414+0-928 
J^-*J^*-0.536 

-1.160 



4.104"' 



2.i 
5.102"' 



1.894"' 



0.718 



-0.75: 
-0.581 



-+0.03 
-0.02 

s+0.06 
-0.05 

-0.59 ±0.06 
-0.59 ±0.04 
-0.48 ±0.07 

-0.86to-^ 



-2.40 ±0.04 
-2.83J 
2.941 
< 



5+0.14 
-0.20 
F0.18 
'-0.25 

-5.0 



o.uy_g 35 

-1 8c;+0.05 
-'-•°>'-0.06 
-3.09 (fixed) 



9 49fi+^-^'^-' 

1+0.697 
-0.658 



3.011"' 



1 Sfi9+0-719 
i.«b^_0.525 



6.439 
3.721 



+ 1.550 
1.230 
1.260 



8.71 



+ CXD 

-4.18 





1016 / 970 


1.62 ±0.03 


979 / 968 


1.60 ±0.03 


976 / 967 




840 / 971 


1.66 ±0.04 


991 / 968 


1 c:4+0.07 
1 Q9+0.20 


886 / 970 
890 / 969 
43 / 118 



CO 
CO 



